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The importance of introducing the SO (5) singlet into the new novel model of non- Fermi liquid 
proposed by LeClair and collaborators is pointed out. This SO(5) singlet order parameter leads to 
a phase transition (or cross-over) between the normal metal and the pseudogap region in the phase 
diagram for the anti-ferromagnetic (AF) phase and the (d-wave) superconducting (SC) phase. A 
non-zero value for this order can increase substantially the critical temperature of the SC transition. 

PACS numbers: 74.25.Dw,74.20Mn,71.27.+a 



It is "well-known that the Landau Fermi liquid model 
does not describe the properties of the high T c super- 
conductivity [H- In search of an alternative model, the 
concept of the renromalization group (RG) flow and Wil- 
son low energy effective field theory should prove useful 
0. Let d be the number of spatial dimensions of the 
system we are interested in. In the action formulation, 
in units where h, ks and the limiting speed are set to 
unity, all terms in the Langrangian density has dimen- 
sion D = d + 1. A contact interaction term (operator) 
with (mass) dimension n < D (so the coupling has di- 
mension D — n) is relevant since it can dominate in the 
low energy limit. Such interactions are relevant opera- 
tors. Interaction operators with dimensions higher than 
D are irrelevant in low energy and so can be ignored in 
leading approximations. Interaction operators with di- 
mensions equal to D are marginal. Their low energy rel- 
evance depends on whether they are asymptotically free 
(becoming strong at low energy scales) or not. 

If we treat the quasi-particle (electron) field \ with 
respect to the Fermi surface as a Dirac-like spin 1/2 
fermion, the kinetic term of \ involves only one deriva- 
tive, so x nas classical dimension of [v] — d/2. The 
4-Fermi interaction has naive classical dimension 2d. In 
d = 3, the effective dimension of the Cooper pair inter- 
action (including the ^-function for momentum conser- 
vation) is reduced from 6 to 4, so the BCS interaction 
is marginal 0]. As is well known, it plays a crucial role 
in low T c (i.e., normal) superconductivity. In the d = 2 
case for cuprates as high T c superconductors, data clearly 
shows that the BCS theory does not apply. The classi- 
cal [x] = d/2 = 1, so the naive 4-Fermi interaction in the 
Landau Fermi liquid model has dimension 4, which makes 
it irrelevant. Apparently we are left with no known rel- 
evant or marginal interaction to include. Such a theory 
leads to free electrons, i.e., normal metal. 

In the novel work by LeClair and collaborators [!, 0, [|| , 
an alternative \ kinetic term with 2 derivatives (of Klein- 
Gordon type) is proposed and analyzed. In this model, 
the [x] = {D - 2)/2 = 1/2 and the 4-Fermi interaction 
(with classical dimension 2) is relevant. Its strong in- 
teraction is crucial in deriving the phase diagram. If 
the model captures the key features of high T c supercon- 
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FIG. 1: The schematic phase diagram as a function of 
x = l/g, which is proportional to the hole doping, 
h(x) (xx — xo- = 1/.9* — 8 is the renormalization 
group fixed point. The y-axis also measures the critical 
temperatures. The gap 5' u (solid dark blue curve) is the 
pseudogap (and/or the spin gap) transition discussed in 
the paper. The dashed line indicates the scaling 
r = A/A c . The presence of 8' u enlarges the gap S' s for 
the anti-ferromagnetic (AF) phase and the gap S' q for 
the d-wave superconducting (SC) phase (solid red 
curves) with respect to that obtained in Ref. 5]. 



ductivity, one can incorporate more structures into it. 
Here our discussion closely follows Ref. [B| , which shows 2 
phase transitions : the anti-ferromagnetic (AF) and the 
(d-wave) superconducting (SC). In this paper, we show 
that there is a natural (automatic) extension of the ba- 
sic concept to include another order parameter, which 
leads to a third phase transition. This is identified with 
the pseudogap transition. In contrast to the AF and SC 
phases, which have spin and charge (respectively) proper- 
ties, the pseudogap phase does not have any distinguish- 
ing feature and may be considered as a cross-over. This 
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completes the construction of the basic model. 

This novel model has a SO(5) symmetry and the 
order parameters in the vector 5 representation has been 
studied in Rcf . [5] . The neutral order parameters in 5 lead 
to the anti-ferromagnetic (AF) phase and the charged 
order parameter in 5 leads to the superconducting (SC) 
phase. In this paper, we include the SO (5) singlet or- 
der parameter. We argue that it is responsible for the 
pseudogap. Its condensate also raises the other critical 
temperatures, in particular that for the superconduct- 
ing transition. In principle, the overall features of the 
phase diagram are essentially parameter-free. (In prac- 
tice, some parameters remain to be determined and oth- 
ers to be better calculated.) Quantitative properties re- 
quires the value of the dimensionless coupling g at the 
energy-momentum cutoff scale A c , <?o — i?(A c ). It is nat- 
ural to take A c to be the inverse of the lattice spacing 
a. To obtain the critical temperatures, one also needs to 
know the Fermi velocity vf- 

Around the Fermi surface at less than half-filling, the 
model is given by 
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dtd 2 x (^2,dtx a d t 
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where a =1,1 and the ± indices on the fields x ± cor ~ 
respond to electric charge. The mode expansions of the 

P ' x), 



fields are (w p = yp 2 and p ■ x = u 

X +(x,t) = J ^f^=(-bi p e-^+a ap e^) 

With {xt-.X^} = ±i5 ap and {Xi,Xj} = 0, the 
(a, b, a',w) are the normal electron and hole annihila- 
tion and creation operators. The Hamiltonian H is 
pseudo-hermitian (H^ = CHC where C is a special 
unitary operator satisfying C^C = C 2 = 1 and distin- 
guishes particles and holes: CaC = a, CbC = —b and 
(x~)^ c = C(x~)^C = x + )i ensuring a unitary time evo- 
lution and real eigenvalues 0, @- Here the relevant 4- 
Fermi contact interaction is unique, up to the strength of 
the coupling g > 0, which has dimension 1. The coupling 
g(A) is scale dependent, so it is convenient to introduce 
the dimensionless coupling 5(A) = g(A)/A = 1/x. Ref. 4] 
finds 
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which has an attractive fixed point at g* = ~ 1/8. 
Suppose g(A c ) = go = 1/xq. Assuming that (3(g) has no 
other fixed points between g$ > g > g*, we have 

A c <?o 



3(A) 



A + 8(A C - A)g 



(2) 



Hole doping goes as h cx x — xq, so we are interested in 
the case where go 3> g* = 1/8 and study the behavior of 
the theory as it flows towards c/» as A — > 0. 

The model has an automatic SO (5) symmetry 0]. 
Treating the xs as the spinor representation 4 of SO (5), 
we see that 4x4x4x4 = 1 + ... where 1 is a SO (5) 
singlet. The interaction operator is such a singlet. Now 
we can write 4x4 = 1 + 5 + 10, where 1, 5 and 10 are 
the singlet, the vector and the adjoint representation of 
SO (5), respectively. That is, one can write the bilinears 
of x fields in terms of pair fields or order parameters. 
The singlet 1 is given by 



and the vector 5 is given by 
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= (x T xl > Xi X T , ^(Xt Xt - Xj. *4 )' X{Xj,Xi X T ) 

where <fi± = (<f> x ± i(f> y )/y/2. In terms of the subgroup 
SU{2) x U{1) of SO(5), the = x _ ^X + /\/2 is an electri- 
cally neutral SU (2) vector and <f)f are Cooper pair fields 
of charge ±2 which are SU (2) spin singlets. Note that 



$ ■ $ = 2<t> + 4>- +4>l- 20 e + C 



(4) 



where the "— " sign follows from pseudo-hermiticity. Due 
to the Grassmanian property of the x s > the interaction 
Hamiltonian can be written in terms of any combination 
of the order parameters, 



//;„, = / d 2 x 8-K 2 g x T x\ X t X; = ' 

d 2 x 8n 2 g(j) 2 z = + J d 2 x 8it 2 g(j)\ 
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However, the interaction Hamiltonian H{ nt can now be 
generalized in a manifestly SO (5) invariant fashion to 
allow more structure to the model : 

d 2 xd 2 x' \$(x) ■ $(x')Vs(x,x') +(j> (x)(t>o(x')Vo(x,x' 

(5) 

where a marginal <fio$ ■ $ term is also possible. Electro- 
magnetic couplings of <jx^ imply that the SO (5) symmetry 
is valid only approximately. Since the x s are Grassma- 
nian, a 10 must be constructed with a derivative, e.g., 
Q ~ XX (they form the SO(5) generators); but its di- 
mension [Q] = 2, so a singlet interaction of the form 
J2QijQij i s irrelevant under naive power counting. (An 
interaction ~ 5 x 5' where 5' contains a derivative will 
be marginal, but such a singlet breaks Lorentz proper- 
ties.) So the model ([5|) is the most general relevant in- 
teraction consistent with the SO (5) symmetry (based on 
naive power counting). However, it is important to note 
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that, under appropriate conditions (as in some momen- 
tum subspaces), higher dimensional operators may be- 
come relevant/marginal and so must be included (e.g., 
Cooper pairs in the BCS model and the d-wave here). 
Some of these higher dimensional SO (5) symmetric op- 
erators (e.g., the d-wave here) would be generated auto- 
matically via quantum corrections. 

Introducing the §\ operator is the new ingredient be- 
yond the analysis of Ref. [|| , so we shall focus on the dy- 
namics of this term. Although classically, both [cf>o] = 1 
and [<fi] = 1, their anomalous dimensions have been eval- 
uated in Ref. [H. At the RG fixed point, [0 O ] ~ 5/8 and 
[(f)] ~ 3/2 in the e-expansion. This suggests that quantum 
effects will make the (j>Q interaction more relevant, while 
pushing the • $ towards marginality. In this sense, the 
4>q interaction is the most relevant and important one. 

Following the standard Hubbard-Stratonovich ap- 
proach, we introduce the auxiliary fields s, q and u 
coupled to the order parameters with the action 

S aux = J dtd 2 x + V2u<fj + \/2s ■ $+ q + 4>~ + q~ <\>t 

1 



1 2 1 ~ ~ 

-u — ——z — s ■ s 



8ir 2 g u 8ir 2 g s 8n 2 g q 
Variations SS aux — imply 



q q 



u = 8ir 2 g u (f> Q /V2, s = 8n 2 gJ/V2, q ± = 8^ g q cf>t , 

Substituting these back into the action reproduces the 
original model if 



g u /2 + 3g s /2 - g q 



(6) 



This relation is valid only classically. SO (5) symmetry 
implies g s — -g q . 

Instead, integrating out the \ fields (in euclidean 
space) yields an effective potential in terms of s, q^ and 
u, 
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where the mass m should be interpreted as the tempera- 
ture m cx T. 

Varying V c g with respect to u, s and q^ yields the gap 
equations respectively, 
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Since u, s = |s| and q = q ± have dimension 2, we in- 
troduce the dimensionless S-j : 



^A 2 , s = S 2 A 2 C and 
q = 5 2 A 2 . Here, the 5j are energy scale (A)-dependent. 
The 3 gap equations are coupled to one another, so a gen- 
eral analysis is a little complicated. The s gap equation 
and the_o ± gap equations have been derived and analyzed 
in Ref.[5j. So let us consider the u gap equation ([7]) here. 

Since quantum effects render the interaction opera- 
tor more relevant than the others, it is reasonable to first 
treat the u gap equation ignoring possible s and q gaps. 
Then one can study the s and q gap equations treating 
the solved u as a background. Let m = s = q = 0, so 
Eq.0 simplifies to 



-4^ 



dp- 



(10) 



Here u — > as g u — > 0; so u — * is the solution as g u — > 0. 
Note that this equation PU|) is different from the other 2 
gap equations in that RHS of Eq. (jTUJ) is not proportional 
to u. For g u > 0, the s-wave attractive interaction yields 
a gap at scale A for S u > 1, 



6 a = ^ ftanh- 1 1/S U A 



(11) 



For large g u /A c , we have 5 U ~ (Agu/SAc) 1 / 4 . For 6 U > 1, 
we have 5 U = 1 + 2e~ Ac / 29u — > 1 if g u decreases as A 
decreases, as expected. Including classical scaling, we 
have S' u ~ 5 u r = S U A/A C . This is shown in Fig 1. 

However, g u may be negative. Following Eq.((6]), one 
may expect g u ~ —g. In this case, we may solve Eq. (|10p 
for u perturbatively. To leading order in g, we now have 
u ~ igA c = 4gAA c In this crude approximation, we have 
5' u = y^u/A 2 ~ y/4gA/A c , where g is approximately 
given in Eq.((2|). Because of the form of the equation 
(fT0|) . this is more like a cross-over than a phase transi- 
tion. Note that irrespective of the sign of g u , we have a 
positive non-zero u solution. 

The presence of a non-zero u does not break the SO (5) 
symmetry, so no Goldstone mode is introduced. Since 
it does not break the SU (2) spin or the electromagnetic 
U(I) symmetries, its detection is more subtle than the 
other 2 phase transitions. It should show up in the prop- 
erties of, e.g., the specific heat and the paramagnetic 
susceptibility [§]. 

The S'j (or Sj) (j = u,s,q) gaps are scale-dependent. 
Analogous to the anomalous mass dimension 7 m for the 
mass, we expect an anomalous gap dimension, jj = 
d Q-^\ , whose 2-loop result 0] may not be that informa- 
tive for strong couplings. Clearly their scaling behaviors 
remains to be determined. So far, except for the run- 
ning of the coupling, classical scaling ( as a function of 
r = A/A c has been used. 

Next, we are interested in the determination of the 
pseudogap temperature Tp, which is related via the u gap 



equation ([7]) by 



0, where the mass m = oT . 
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Q estimates that a ~ 1.7. This leads to the critical 
temperature for the S'j gap, 

Tj ^ 9(5} j = u,s,q (12) 

where Cj is some normalization factor. In Fig. 1, we see 
that all of the AF and SC regions are probably inside 
the region with a non-zero u value, so we expect those 
regions will be modified by the presence of u. 

The s gap equation (jHJ) is attractive for weak doping, 
so it has a solution, with a non-zero 5' s that breaks the 
SU(2) symmetry [f|. It is shown in Fig. 1. 

With g q positive, the s-wave g r± gap equation © has a 
repulsive potential and so does not have a solution. Phys- 
ically, it is not surprising since the pairing of 2 electrons 
has a repulsive Coulomb force between them, which is 
absent in the u and s pairings. This follows from Eq.([S]). 
which comes from the "— " sign in Eq. ([?]), a consequence 
of the definition of cj)^, which in turn is dictated by 
pseudo-hermiticity. Introducing non-constant auxiliary 
pair fields and one-loop contribution, Kapit and LeClair 
derived a momentum-dependent gap equation, which has 
a d-wave attractive channel for Cooper pairs This 
yields the superconducting phase, i.e., the SC region in 
the phase diagram schematically shown in Fig. 1. The ef- 
fective d-wave coupling gi = 4<? 2 /25A 3 decreases rapidly 
as doping is decreased (i.e., A increases). Once gi be- 
comes too small, the potential is too weak to bind, so the 
SC region terminates. This is the picture in the absence 
of a u value. Now we like to make a crude estimate of 
the impact of a non-zero u on the SC region. 

The high doping end of the SC dome region is termi- 
nated at the critical fixed point x* = 8. The upper SC 
transition at doping h(x* = 8) = 3/2n 2 ~ 0.15 is second- 
order and corresponds to the fixed point of the renor- 
malization group. This is universal in this model. So we 
expect the main effect of a positive u is to raise the height 
of the SC dome, that is, raising the critical temperature 
T$c for the superconducting transition. Since the lower 
doping end is due to the weakness of the d-wave coupling, 
the low doping end (x\) may decrease some, but we do 
not expect it to change much. 

The (maximum) critical temperature Tn = Tsc( u = 0) 
in the absence of u is estimated in Ref.[5j, in which the 
reader can find the details. Restoring in Eq. (fT2|) h and ks 
and the Fermi velocity vp ~ 210 km/s as the limiting 
speed, and putting in values obtained in Ref. c q ~ 0.5, 
a ~ 1.7, S' qQ — 6' g (u = 0) ~ 0.11 and A c — 1/a, where the 
lattice spacing is a ~ 3.8 x 10~ 10 to, one obtains 

Tsc (u = 0) = , g ° = Jc q 6' q0 ~ U0°K (13) 

a kb 

where the temperature J = vph/aaks — 2500°X. The 
actual Neel temperature is lower than J, presumably due 
to inter layer couplings. 



Looking at the effective mass term in the denominator 
of the q gap equation (JOj) , we see that m 2 — > m 2 — u. 
This means that, in the presence of a non-zero 5' u , the 
q gap equation can go to a higher T without losing its 
gap solution. That is, the critical temperature Tsc is 
enhanced. The effect is largest for largest A (smallest 
doping) where gi has not yet cut off the q gap solution. 
Taking go — * 00 and S' u ~ \Jr ~ 0.2, 

T SC = J [(c q S' q0 ) 2 + (c u 6' u ) 2 } 1/2 ~ 280°K (14) 

where we simply take c u ~ c q . Although this value of 
Tsc is for illustrative purpose only, we see that a u > 
value (irrespective of the sign of g u ) enhances the Tsc- 
At this doping value, the pseudogap transition temper- 
ature Tp = Jc u S' u ~ 250°K is bigger than To. Since 5' u 
gap vanishes above Tp, we expect Tsc to loosely track Tp 
until it reaches the lower doping value where the critical 
temperature drops to zero (i.e., the q gap disappears). As 
pointed out in Ref.Q, decreasing the lattice spacing or 
increasing the Fermi velocity increases the critical tem- 
peratures. We see that a more accurate determination of 
the scaling properties of the S'j (as well as a and Cj ) will 
provide a crucial test of the model. 

Similarly, the combination of q + q~ — s 2 in the denomi- 
nators inside the gap equations (JH [9|) suggests that bind- 
ing in one channel strengthens the binding in another 
channel (this may happen if g(A c ) < g*, or in the elec- 
tron doped region). Instead of competing, the order pa- 
rameters tend to strengthen each other. In summary, the 
inclusion of the SO (5) invariant order completes the con- 
struction of the basic model. The interplay of the 3 order 
parameters yields a non-trivial phase diagram. 
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